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WHAT IS CLAIMED IS: 

Please amend the claims as follows: 

1 . (Currently Amended) A numerical calculation method for a physical 
quantity U practiced on a computer by solving A-U = f 5 wherein A is a coefficient 
matrix (in N rows by N columns; wherein N is a positive integer) obtained 
through discrete of a partial differential equation to be satisfied by said physical 
quantity U, and f is an inhomogeneous term (source term), comprising the 
processes of: 

setting an initial value U° of said physical quantity U; 

setting 0 as an initial value of a number m of repeating times, giving 0 as 
an initial value of a perturbation quantity 0and setting (f - A.U°) as an initial 
value r° of a residual r; and 

repeatedly executing a first step and a second step while incrementing 
said number m of repeating times until an approximate solution U m is converged, 
wherein said first step includes the steps of: 

obtaining a predicted approximate value i|/ m of A* <j> = r m through 
repeated calculation performed by a first calculation unit including an internal 
solver, and 

said second step includes the steps of: 

obtaining, from said predicted approximate value \|/ m , a corrected 
approximate value (j> m for minimizing L 2 norm of a residual r^ r" 1 * 1 t hrough an 
optimization routine performed by a second calculation unit; and 

giving (U m + <f> m ) as an approximate solution U m+1 and giving (r m - 
A-<zT) as a residual r^ 1 , 

wherein in said second step, obtained elements of a vector sequence 
A • <p m are sampled by a given sampling method to be stored in a memory, and 

a residual minimization coefficient ai m (wherein 1 = 1, . . ., L) used for 
obtaining said corrected approximate value <p m is approximately obtained by 



NVA298719.1 



- 5 - Docket No.: 740819-1036 

using elements of a vector sequence A - <j> v (wherein k = m - L + 1, . . ., m - 1) 

stored in said memory. 

2. (Original) The numerical calculation method of Claim 1, 
wherein in sampling of elements bi, b 2 , - . . and b N of said vector 

sequence A- <p m performed in said second step, elements b\ (wherein i efl) are 

selected, whereas a subset £1 is defined as follows: 
ft = {i : mod [i, lg] = 1} U {i : |fi/a u | > P} 
wherein lg is an integer, P is a real number, £ is an element of said 

source term and a^ is a diagonal term on the ith row in the ith column of said 

matrix A. 

3. (Currently Amended) A numerical calculator for a physical quantity U 
by solving A-U = f, wherein A is a coefficient matrix (in N rows by N columns; 
wherein N is a positive integer) obtained through discrete of a partial differential 
equation to be satisfied by said physical quantity U, and f is an inhomogeneous 
term (source term), performing the processes of: 

setting an initial value U° of said physical quantity U; 

setting 0 as an initial value of a number m of repeating times, giving 0 as 
an initial value of a perturbation quantity ^and setting (f - A.U°) as an initial 
value r° of a residual r; and 

repeatedly executing a first step and a second step while incrementing 
said number m of repeating times until an approximate solution U m is converged, 
wherein said first step includes the steps of: 

obtaining a predicted approximate value V|/ m of A- <f> = r m through 
repeated calculation performed by a first calculation unit including an internal 
solver, and 

said second step includes the steps of: 

obtaining, from said predicted approximate value \|/ m , a corrected 
approximate value <f> m for minimizing L 2 norm of a residual r"*r m+1 through an 
optimization routine performed by a second calculation unit; and 

giving (U m + ^ m ) as an approximate solution u m+1 and giving (r m - 
A-0 m ) as a residual r m+1 , 
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wherein in said second step, obtained elements of a vector sequence 
A- <f> m are sampled by a given sampling method to be stored in a memory, and 

a residual minimization coefficient a™ (wherein 1 = 1, L) used for 
obtaining said corrected approximate value <j> m is approximately obtained by 
using elements of a vector sequence A- <f> k (wherein k = m - L + 1, . . ., m - 1) 
stored in said memory. 

4. (Original) The numerical calculator of Claim 3, 

wherein in sampling of elements bi, b 2 , ... and b N of said vector 
sequence A- <j) m performed in said second step, elements bi (wherein i gQ) are 
selected, whereas a subset Q. is defined as follows: 
Q = {i : mod [i, lg] = 1} U {i : W*A > P> 

wherein lg is an integer, p is a real number, £ is an element of said 
source term and aji is a diagonal term on the ith row in the ith column of said 
matrix A. 

5. (Currently Amended) A recording medium that stores a numerical 
calculation program for a physical quantity U by allowing a computer to solve 
A U = f, wherein A is a coefficient matrix (in N rows by N columns; wherein N 
is a positive integer) obtained through discrete of a partial differential equation to 
be satisfied by said physical quantity U, and f is an inhomogeneous term (source 
term), 

wherein said numerical calculation program makes said computer to 
execute the processes of: 

setting an initial value U° of said physical quantity U; 

setting 0 as an initial value of a number m of repeating times, giving 0 as 
an initial value of a perturbation quantity <f> and setting (f — A.U°) as an initial 
value r° of a residual r; and 

repeatedly executing a first step and a second step while incrementing 
said number m of repeating times until an approximate solution U m is converged, 
wherein said first step includes the steps of: 
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obtaining a predicted approximate value \|/ m of A- <j> = r m through 
repeated calculation performed by a first calculation unit including an internal 
solver, and 

said second step includes the steps of: 

obtaining, from said predicted approximate value \|/ m , a corrected 
approximate value ^ m for minimizing L 2 norm of a residual f"* r m+1 through an 
optimization routine performed by a second calculation unit; and 

giving (U m + (f> m ) as an approximate solution U m+1 and giving (r m - 

A- <*T) as a residual r m+ \ 

wherein in said second step, obtained elements of a vector sequence 
A- ^ m are sampled by a given sampling method to be stored in a memory, and 

a residual minimization coefficient a™ (wherein 1 = 1, . . ., L) used for 
obtaining said corrected approximate value <f> m is approximately obtained by 
using elements of a vector sequence A- <f> v (wherein k = m - L + 1, . . m - 1 ) 
stored in said memory. 

6. (Original) The recording medium of Claim 5, 

wherein in sampling of elements bi, b 2 , ... and bu of said vector 
sequence A- <j> m performed in said second step, elements bi (wherein i efi) are 
selected, whereas a subset Q is defined as follows: 
Q = {i : mod [i, l g ] = 1 } U {i : \ffa\ > p} 

wherein lg is an integer, p is a real number, £ is an element of said 
source term and a** is a diagonal term on the ith row in the ith column of said 
matrix A. 



NVA298719.1 



